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The concept of gauge invariance can be considered one of the most subtle and useful concept in 
theoretical physics since it can permit the comprehension of difficult systems in physics with an 
arbitrary choice of a reference frame at every instant of time. It is always desirable to have a bridge 
between gauge invariant and noninvariant theories. Once established, this kind of mapping between 
first-class (gauge invariant) and second-class systems, in Dirac's formalism can be considered as a 
sort of duality. In this paper we investigate this "duality" obtaining a gauge invariant theory starting 
with a noninvariant one. We analyzed both Abelian and non-Abelian theories and the procedure 
used is the recent dual (also called symplectic) embedding formalism. We believe that this method 
is the most convenient one since it is not plagued by the ambiguity problems that torments BFFT 
and other iterative methods. We demonstrated exactly that this "dualization" method used here 
does not require any special modification to handle with non-Abelian systems, which is also a new 
result described in this paper. To prove the gauge invariance we used just the Dirac constraint 
technique. It is relevant to say here that, although this work is lengthy, the majority of the results 
presented here are new in the literature. The results that are not new were reproduced within a 
new perspective. 
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which can be thought as one in which the dynamical vari- 
ables are determined with relation to a reference frame 
whose choice is arbitrary at every instant of time. The 
relevance of a gauge theory, in few words, is that the 
physically important variables are those that are inde- 
pendent of the local reference frame [2;]. Whenever a 
change in the arbitrary reference frame causes a trans- 
formation of the variables involved we have the so-called 
gauge transformation. This physical variables are then 
well known as gauge invariant variables. Such gauge the- 
ories and gauge transformations are the cornerstones of 
the construction of the standard model, to mention only 
one of its successful applications. We will talk with more 
detail about gauge theories and quantization in a mo- 
ment. 

Since the Hamiltonian formulation is considered by 
many as the more fundamental formulation of a physical 
theory we use the Dirac brackets in order to show exactly 
the gauge invariance of the final actions obtained. 

Second-class — ^ first-class. Gauge theories have 
played an important role in field theories since they are 
related with the fundamental physical interactions of Na- 
ture. In a more general sense, those theories have gauge 



I. INTRODUCTION 

This work deals with the transformation of second- 
class systems into gauge theories. Although there is an 
extensive literature about this subject, we obtain some 
results here that are altogether new. In a self-sustained 
and pedagogical way we analyze gauge invariance and 
dual embedding issues. 

A map with these conversion features can be under- 
stood as a kind of duality since we will show precisely that 
the final action have only first-class constraints. This du- 
ality connection characterizes both systems as represent- 
ing the same physical properties. We will discuss these 
ideas carefully later. 

As well known, after the procedure, having only first- 
class constraints, the final system is a gauge theory. 
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symmetries defined by some relations called, in Dirac's 
language, first-class constraints The quantization of 
these theories demands a special care because the pres- 
ence of gauge symmetries indicate that exist some su- 
perfluous degrees of freedom, which must be eliminated 
(before or after) with the implementation of a valid quan- 
tization process. 

The quantization of first-class systems was formulated 
both in Dirac's [l| and path integral Q point of view. 
Later on, the path integral analysis was extended by 
Batalin, Fradkin and Vilkovisky [4] in order to preserve 
the BRST symmetry [|. 

On the other hand, the covariant quantization of 
second-class systems is, in general, a difficult task be- 
cause the Poisson brackets are replaced by Dirac brack- 
ets. At the quantum level, the variables become opera- 
tors and the Dirac brackets become commutators. Due 
to this, the canonical quantization process is contami- 
nated with serious problems such as ordering operator 
problems and anomalies [7] in the context of non- 
linear constrained systems and chiral gauge theories, re- 
spectively. In view of this fact, it seems that it is more 
natural and safe to work out the quantization of second- 
class systems without invoking Dirac brackets. Actually, 
it was the strategy followed by many authors over the 
last decades. The noninvariant system has been embed- 
ded in an extended phase space in order to change the 
second-class nature of constraints into first-class. 

In this way, the entire machinery Q for quantiz- 
ing first class systems can be used. To implement this 
concept, Faddeev suggests to enlarge the phase space 
with the introduction of new variables to linearize the 
system, which were named after, as the Wess-Zumino 
(WZ) variables. This idea has been embraced by many 
authors and some methods were proposed and some 
constraint conversion formalisms, based on the Dirac 
method [T|, were constructed. Among them, the Batalin- 
Fradkin-Fradikina-Tyutin (BFFT) [10] and the iterative 
[Tlj methods were strong enough to be successfully ap- 
plied to a great number of important physical systems. 
Although these techniques share the same conceptual ba- 
sis iOl and follow the Dirac framework [1], these constraint 
conversion methods were implemented following differ- 
ent directions. Historically, both BFFT and the iterative 
methods were introduced to deal with linear systems such 
as chiral gauge theories [111, [13 in order to eliminate the 
gauge anomaly that hampers the quantization process. 

In spite of the great success achieved by these meth- 
ods, they have an ambiguity problem [T3| . This problem 
naturally arise when the second-class constraints are con- 
verted into first class ones with the introduction of WZ 
variables. Due to this, the process of constraint conver- 
sion process may become a hard task, as shown in (Tsj . 
Duality. Duality is a very useful concept in field theory 
and statistical mechanics since there are very few analytic 
tools available for studying non-perturbative properties 
of systems with many degrees of freedom. Always stud- 
ied as one of the first applications of duality basics, the 



electromagnetic duality is a pedagogical and interesting 
example for the interested reader 14j. 

Recently, the so-called gauging iterative Noethcr dual- 
ization method 15] has been shown to thrive in establish- 
ing some dualities between models [16 - 18] . This method 
is hinges on the traditional concept of a local lifting of a 
global symmetry and may be realized by an iterative em- 
bedding of Noether counterterms. However, this method 
provides a strong suggestion of duality since it has been 
shown to give the expected result in the paradigmatic du- 
ality between the so-called self-dual model [11] and the 
Maxwell-Chern-Simons theory in three dimensions dual- 
ity. This correspondence was first established by Deser 
and Jackiw (20| and using the parent action approach 

The paper. We have organized this paper as follows. 
In section 2, we review the dual embedding formalism 
in order to settle the notation and familiarize the reader 
with the fundamentals of the formalism. In section 3, 
we will begin to make the application of the ideas dis- 
cussed before in some Abelian models. We initiate with 
the Proca model in order to set up the general ideas dis- 
cussed in Section 2. We will see that a different zero- 
mode can bring us an equivalent action different from 
the literature. After that, we apply this formalism in 
two important models in high energy physics. The first 
one, in section 4, is the nonlinear sigma model (NLSM) 
[2^ , which is an important theoretical laboratory to learn 
the basics about asymptotically free field theories, as 
dynamical mass generation, confinement, and topologi- 
cal excitations, which is expected in the realistic world 
of four-dimensional non-Abelian gauge theories. The 
second, in section 5, is the bosonized chiral Schwinger 
model (CSM), which has attracted much attention over 
the last decade, mainly in the context of string theo- 
ries [131, also due to the huge progress in under- 
standing the physical meaning of anomalies in quantum 
field theories achieved through the intensively study of 
this model. Through this section, we will first compute 
the Dirac brackets among the phase space fields and af- 
ter that, the gauge invariant version of the model will 
be obtained. Later on, the gauge symmetry will be in- 
vestigated from the Dirac point of view. Section 6 is 
devoted to an application of the formalism to the non- 
Abelian Proca model. In this section, we will show that 
this gauge-invariant formalism does not require modifi- 
cations to deal with non-Abelian models as demanded by 
the BFFT method. In section 7, we analyze the symme- 
tries of the rotational fluid model with a new extra term, 
which introduce a dissipative force into the system. The 
objective is to promote an approximation to reality where 
we always have dissipation. Our concluding observations 
and flnal comments are given in Section 8. 
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II. THE DUAL EMBEDDING FORMALISM 

The formalism is developed based on the symplectic 
framework [2^, , that is a modern way to handle with 
constrained systems. The basic object behind this for- 
malism is the symplectic matrix, i. e., if this matrix is 
singular, the model presents a symmetry, if it is not sin- 
gular, the Dirac brackets can be obtained. In this way, 
we propose to render nonsingular symplectic matrix to a 
singular one. It will be carried out introducing arbitrary 
functions that depend on the original and WZ variables 
into the first-order Lagrangian. To appreciate this point, 
a brief review of the symplectic formalism will be fur- 
nished and, after that, general ideas of the symplectic 
gauge-invariant formalism will be presented. This for- 
malism, differently from the BFFT and other iterative 
constraint conversion methods, does not require special 
modifications into the procedures to convert Abelian or 
non-Abelian set of second-class constraints into first-class 
ones. 

This technique follows the Faddeev-Shatashivilli's sug- 
gestion |9|] and is set up on a contemporary framework to 
handle constrained models, namely, the symplectic for- 
malism [2^. In the following lines, as we said above, we 
will try to keep this paper self-sustained reviewing the 
main steps of the dual embedding formalism. We will 
follow closely the ideas contained in [251] . 

Let us consider a general noninvariant mechanical 
model whose dynamics is governed by a Lagrangian 
C{ai,di,t), (with i = 1,2,...,A^), where Oi and di are 
the space and velocity variables, respectively. Notice 
that this model does not result in a loss of generality nor 
physical content. Following the symplectic method the 
zeroth-iterative first-order Lagrangian one-form is writ- 
ten as 



and the symplectic variables are 



;(0)e 



Gi, with 61= 1,2,..., 

Pi, with 6* = A-t- 1, A-l- 2, . 



,,2A, 



(2.1) 



(2.2) 



where Ag are the canonical momenta and is the 
symplectic potential. From the Euler-Lagrange equations 
of motion, the symplectic tensor is obtained as 



dA, 



(0) 



If the two-form 



/ 



dA^ 



^fepd^ A d^^ 



(2.3) 



is singular, the symplectic matrix (|2.3I) has a zero-mode 
(i/*^°)) that generates a new constraint when contracted 
with the gradient of the symplectic potential. 



This constraint is introduced into the zeroth-iterative 
Lagrangian one-form equation (j2.ip through a Lagrange 
multiplier ry, generating the next one 

C^'Ut = Af^d^^"^' + dryri^o) - V'^''\Odt, 

= _ v^(i)(^)rfi^ (2.5) 

with 7 = 1,2,..., (2A + 1) and 



4i) = (4°\l7(o)). 



(2.6) 



As a consequence, the first-iterative symplectic tensor is 
computed as 



dA. 



(1) 



5^(1)7 9^(l)/3 



(2.7) 



If this tensor is nonsingular, the iterative process stops 
and the Dirac brackets among the phase space variables 
are obtained from the inverse matrix (/^Jj')"^ and, con- 
sequently, the Hamiltonian equation of motion can be 
formulated and solved, as discussed in It is well 

known that a physical system can be described at least 
classically in terms of a symplectic manifold A^. From 
a physical point of view, Ai is the phase space of the 
system while a nondegenerate closed 2-form / can be 
identified as being the Poisson bracket. The dynamics 
of the system is determined just specifying a real-valued 
function (Hamiltonian) H on the phase space, i.e., one of 
these real-valued function solves the Hamiltonian equa- 
tion, namely. 



t(A)/ = dH, 



(2.8) 



and the classical dynamical trajectories of the system in 
the phase space are obtained. It is important to mention 
that if / is nondegenerate, the equation (|2.8p has a very 
unique solution. The nondegeneracy of / means that the 
linear map b : TM T*M defined by b(A) := \>{X)f is 
an isomorphism. Due to this, the equation (|2.8p is solved 
uniquely for any Hamiltonian (A = [>~^{dH)). On the 
other hand, the tensor has a zero-mode and a new con- 
straint arises, indicating that the iterative process goes 
on until the symplectic matrix becomes nonsingular or 
singular. If this matrix is nonsingular, the Dirac brack- 
ets will be determined naturally. 

In [2^ , the authors consider in detail the case when / is 
degenerate. The main idea of this embedding formalism 
is to introduce extra fields into the model in order to 
obstruct the solutions of the Hamiltonian equations of 
motion. We introduce two arbitrary functions that hinge 
on the original phase space and on the WZ variables, 
namely, ^/{ai,pi) and G{ai,pi, rf), into the first-order one- 
form Lagrangian as follows 



4°'d?'"^'' + *d?7-V''"^(Orf^, (2-9) 



(2.4) 
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with 



(2.10) 



where the arbitrary function G{ai,pi,r]) is expressed as 
an expansion in terms of the WZ field, given by 

oo 

G(a„p„77) = ^g(")(a„p„77),, (2.11) 



n=l 



where 



and satisfies the foUowing boundary condition 

G(a„p„77 = 0) = 0. (2.12) 

The symplectic variables were extended to also encom- 
pass the WZ variable ^(°)^" = (C(°)^77) (with 9 = 
1,2,..., 2N + 1) and the first- iterative symplectic poten- 
tial becomes 

oo 

V^°^ {ar,Pr, v) = V("'> {a,,p,) + J2 5'"^ V)- (2-13) 

n=l 

In this context, the new canonical momenta are 



^(0) ^ I A^g°\ with 6 =1,2,. . . ,2N 
s I with 6i= 2N + 1 

and the new symplectic tensor is given by 



(2.14) 



that is 



'0/3 Q^(0)6 Q^(Q)I3^ 



/AO) .(0)N 
f(0) _ I ■1013 J0ri 



(2.15) 



(2.16) 



To sum up, basically, we have two steps: the first one 
is addressed to compute \I'(ai,pi) while the second one 



is dedicated to the calculation of G{ai,pi,r]). In order 
to begin with the first step, we impose that this new 
symplectic tensor (/*•"-') has a zero-mode i>, consequently, 
we obtain the following condition 



,io)em = 



(2.17) 



At this point, / becomes degenerated and in consequence, 
we introduce an obstruction to solve the Hamiltonian 
equation of motion given by equation ()2.8|) . Assuming 
that the zero-mode i/^")^ is 



z>(0) = (/ 1) , 



(2.18) 



and using the relation given in (j2.17p together with 
(|2.16p . we have a system of equations, 



where 



„e AO) AO) _ Q 



„(o) _ 



(2.19) 
(2.20) 



The matrix elements ^ are chosen in order to disclose 
the desired gauge symmetry. Note that in this formal- 
ism the zero-mode the gauge symmetry genera- 
tor. It is worth to mention that this feature is impor- 
tant because it opens up the possibility to disclose the 
desired hidden gauge symmetry from the noninvariant 
model. From relation (12.17^ some differential equations 
involving '^{ai^pi) are obtained, i. e., the equation p.l9p . 
and after a straightforward computation, "^{oi^pi) can be 
determined. 

In order to compute G{ai,pi,rj) following the second 
step of the method, it is mandatory that no constraints 
arise from the contraction of the zero- mode (i>('^)^) with 
the gradient of the potential y'^°^(ai,pi, 77). This condi- 
tion generates a general differential equation, which reads 
as 



0, 



9e(i)(a„K,7;) , dg(^\a,,p,,ii) 



g^io)e 



Q^{0)6 



dr] 



drj 



(2.21) 

, 
(2.22) 



that allows us to compute all the correction terms 
Q'^"'\ai^Pi,'q) as functions of 77. Notice that this poly- 
nomial expansion in terms of 77 is equal to zero. Conse- 



quently, all the coefficients for each order in 77 must be 
identically null. Given this, each correction term as func- 
tion of rj is determined. For a linear correction term, we 
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have 



where the symplectic potential is 



For a quadratic correction term, we have 



0. 



(2.23) 



(2.24) 



From these equations, a recursive equation for n > 2 can 
be chosen as, 



9^(0)6 



dr] 



(2.25) 



which permits us to compute the remaining correction 
terms as functions of 77. This iterative process is suc- 
cessively repeated until p.2ip becomes identically zero. 
Consequently, the extra term G{ai,pi,r]) is obtained ex- 
plicitly. Then, the gauge invariant Hamiltonian, identi- 
fied as being the symplectic potential, is obtained from 



(2.26) 



and the zero-mode £>('')^ is identified as the generator of 
an infinitesimal gauge transformation, given by 



(2.27) 



where e is an infinitesimal parameter. 

In the following sections, we will apply the symplec- 
tic gauge-invariant formalism in some second-class con- 
strained Hamiltonian systems, Abelian and non-Abelian. 
We will see that the new results of the following sections 
demonstrate exactly that the dual (symplectic) embed- 
ding method can be seen as a mapping between the orig- 
inal action the final one. 



III. THE ABELIAN PROCA MODEL 

The analysis of the following model clarify the physics 
and gives a deeper insight into the general formalism de- 
scribed in the last section. To this end, let us start with 
a simple Abelian case which is the Proca model whose 
dynamics is ruled by the Lagrangian density. 



1/(0) 



with 



2-^o(5,^^ + ^mMo), 
(3.3) 



TT,: = Ai - diA, 



1^0 



where di = and the dot denote space and time 



derivatives, respectively. The symplectic coordinates 
are ^i*^^ — [Ai^'Ki^ Aq) with the corresponding one-form 
canonical momenta given by 



(0) 

a\l = 7r,j, 



,(0) 



JO) 



(3.4) 



The zeroth-iteration symplectic matrix is 



p(0) _ 









I 6'^^\x-y), 




(3.5) 



which is a singular matrix. It has a zero-mode that gen- 
erates the following constraint. 



17 = a^TT* -I- m^Ao, 



(3.6) 



identified as being the Gauss law. We will include this 
constraint into the canonical part of the first-order La- 
grangian C^^'^ in p.2p introducing a Lagrangian multi- 
plier (/3). The first- iterated Lagrangian can be written 
in terms of £,a^ = ( , TTj , , /3) as 

£(1) =7rM, + 17^-1^(1), (3.7) 
with the following symplectic potential. 



= ^t:^ + + i {Ao' + A,') Ao^. (3., 
The first-iterated symplectic matrix, computed as 



£ = - -F^.F'^'' + - m2 A^^A^, (3.1) 

where m is the mass of the A^ field, g^i, = diag{-\ ) 

and Ffj,i, = d^A^ — d,jA^. Observe that, as well known, 
the mass term breaks the gauge invariance of the usual 
Maxwell's theory. Hence, the Lagrangian density above 
represents a second-class system. 

To perform the dual embedding formalism the La- 
grangian density is reduced to its first-order form as 



1 ° 




























2 


V 




— 


/ 



5^^\x-y), (3.9) 



is a nonsingular matrix. Consequently, the Proca model 
is not a gauge invariant field theory. As settle by the 
method, the Dirac brackets among the phase space fields 
are obtained from the inverse of the symplectic matrix, 
namely. 



£(0) =^M, (3.2) 
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with 



{Mx),A,{y)r ^0, 
{TTi{x),Trj{y)}* ^ 0, 



(3.10) 
(3.11) 



with the fohowing Hamihonian, 



ii?2. + -y-n^A^ 
^^13 + 2 ' 



iTT.MjTT^" + If^. + ^^m'Al (3.12) 



where the phase space metric is 



(3.13) 



which completes the noninvariant analysis. 

At this point we are ready to carry out the symplectic 
gauge-invariant formulation of the Abelian Proca model 
in order to disclose the gauge symmetry present in the 
model. To this end, we will extend the symplectic gauge- 
invariant formalism [2^ , proposed by three of us in order 
to unveil, at that time, the gauge symmetry present on 
the Skyrme model. The basic concept behind the ex- 
tended symplectic gauge-invariant formalism lives on the 
extension of the original phase space with the introduc- 
tion of two arbitrary functions, ^ and G, where both 
rely on both the original phase space variables and the 
WZ variable {6). The former (^f) is introduced into the 
kinetic sector and, the later (G), within the potential sec- 
tor of the first-order Lagrangian. The process starts with 
the computation of 5* and finishes with the calculation 
of G. 

In order to reformulate the Proca model as a gauge 
invariant field theory, we will start with the first-order 
Lagrangian £(0), given in Eq. p.2p . with the arbitrary 
terms, given by. 



£(0) =^«i, + ^*-y(o), 



(3.14) 



with 



1 



1 



1 



1 



- -7T,^ + -Fl^ + -m^A,^-Ao{d,7T' + ^m^Ao) + G, 

(3.15) 

where * ee 'i'{Ai,TTi,Ao,9) and G = G{A,,Tri, Ao,d) are 
the arbitrary functions to be determined. Now, the sym- 
plectic fields are ^i"'' = {Ai,TTi, Aq,9) while the symplec- 
tic matrix is 



f(0) 











Sji 

















a*x 


_a*x 











0*. 
dA 

a* 

dA: 



5(3)(f_y), (3.16) 



39^ d9y 



(3.17) 



where 9^ = e{x), Oy = 9{y), = ^{x) and = ^{y). 

In order to unveil the hidden U{1) gauge symmetry in- 
side the Proca model, the symplectic matrix above must 
be singular, then, vj/ = {Ai,TTi,9). As established by the 
symplectic gauge-invariant formalism, the corresponding 
zero-mode iy^^\x), identified as being the generator of 
the symmetry, satisfies the following relation. 



d'yi^^^\x) f^p{x-y)^Q, 



(3.18) 



producing a set of equations that allows to determine ^' 
explicitly. At this point, it is very important to notice 
that the extended symplectic gauge-invariant formalism 
opens up the possibility to disclose the gauge symmetry 
of the physical model. The zero-mode does not generate 
a new constraint, however, it determines the arbitrary 
function 5* and consequently, obtain the gauge invariant 
reformulation of the model. We will scrutinize the gauge 
symmetry related to the following zero-mode. 



(a, 



(3.19) 



Since this zero-mode and the symplectic matrix p.l6p 
must satisfy the gauge symmetry condition given in Eq. 
p.lSp . a set of equations is obtained and after an inte- 
gration, 5* is computed as 



Hence, the symplectic matrix becomes 



(3.20) 



f(0) 





5,., -af I 

I 

9J 



y), (3.21) 



which is singular by construction. Due to this, the first- 
order Lagrangian is 



(3.22) 



with given in Eq. (I3.15p . 

Now, we start with the second step of the formalism to 
transform the model into a gauge theory. The zero-mode 
v^'^^ does not produce a constraint when contracted with 
the gradient of the symplectic potential, namely. 



.(0) 



(3.23) 
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on the contrary, it produces a general equation that al- 
lows us to compute the correction terms in 9 enclosed 
into G{Ai,TTi, Ao,6), given in Eq. ()2.21|) . To compute 



the correction term linear in 0, namely, Q^-^\ we pick up 
the following terms from the general relation (|2.2ip given 

by, 



( m^A\x)S^^^x-w) 



dF'^{x) \ , dQ^^\x) 
dA\w) 



de{w) 



= 0. 



(3.24) 



r 



After a straightforward calculation, the correction term 
linear in is, 



m^d'AiO. 



(3.25) 



Substituting this result into the symplectic potential 
p.lSp . we have that 

= + \f!^ + i m2 A,^ - A,{d,^^ + \m^ A,) 



- m'd'Ai 



(3.26) 



However, the invariant formulation of the Proca model 
was not obtained yet because the contraction of the zero- 
mode (|3.19p with the symplectic potential above does not 
generate a null value. Due to this, higher order correction 
terms in 9 must be computed. For the quadratic term, 
we have. 



dg^^\x) 

d9{w) 

and after a direct calculation, we can write that. 



(~n?9[x)di5^'^\x-w, 



T(2) 



1 



{d,9f 



= 0, 
(3.27) 

(3.28) 



Then, the first-order Lagrangian becomes. 



of the infinitesimal gauge transformation {50 = ei/*^"^) 
given by. 



= - 

SlTi = 0, 

SAo = 0, 
69 = e. 



die, 



(3.31) 



where e is an infinitesimal time-dependent parameter. In- 
deed, for the above transformations the invariant Hamil- 
tonian, identified as being the symplectic potential V'^^\ 
changes as 



(3.32) 



Now we will investigate the result from the Dirac point 
of view. The chains of primary constraints computed 
from the Lagrangian p.29p are 



Xi 



TTg. 



(3.33) 



Next, these constraints will be introduced into the invari- 
ant Hamiltonian p.30p through the Lagrange multipliers 
and so it can be rewritten as 



jC = Tr' A, + 9^ -V'^"\ 
where the symplectic potential is 



(3.29) 



(0) 



primary 



+ Ai(/)i -H7iXi- 



(3.34) 



The time stability condition for the primary constraint 
requires a secondary constraint, such as 



1 9 

4 »j ■ 2 
1 



2 /I 2 

- m Ai 



m^d'Ai9 + -m^ {d^9f . 



^o) 
(3.30) 



Since the contraction of the zero- mode (p(°)) with the 
symplectic potential above does not produce a new con- 
straint, a hidden symmetry is revealed. 

To complete the gauge invariant reformulation of the 
Abelian Proca model, the infinitesimal gauge transfor- 
mation will be computed also. In agreement with the 
symplectic formalism, the zero-mode D^°^ is the generator 



02 = 9i7r' -I- m'^Ao, 



(3.35) 



and no more constraint appears from the time evolution 
of 02. Now, the total Hamiltonian is written as 



(0) 



total 



V"^°^ -I- Ai0i -I- A202 + 7iXi- 



(3.36) 



Since the time evolution of 0i just allows us to obtain 
the Lagrange multiplier A2, and the constraint xi has 
no time evolution xi = 0, no more constraints arise. 
Hence, the gauge invariant model has three constraints 
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((/>!, (/)2, Xi)- The nonvanishing Poisson brackets among 
these constraints are, 



satisfying a nonhnear constraint tr^ = 1. This model has 
its dynamics governed by the Lagrangian density 



{0i(a;),02(y)} = -m2<5(3)(£_y). 
The Dirac matrix given by 



(3.37) 



C 




(3.38) 



is singular, indicating that the model has indeed, a gauge 
symmetry. However, it also has some nonvanishing Pois- 
son brackets among the constraints, suggesting that the 
model has both first and second-class constraints. It is 
easy to check that xi is a first class constraint and (/>i 
and (j)2 are second-class constraints. In accordance with 
the Dirac method, the set of second-class constraint must 
be taken equal to zero in a strong way, generating then 
the primary Dirac brackets among the phase space fields, 
given by 

{e{x),7rg{y)r ^S^'\x~y). (3.39) 

The gauge invariant version of the Abelian Proca model 
is then governed by the following invariant Hamiltonian, 



(3.40) 



where /i = 0, 1 and "a" is an index related to the 0{N) 
symmetry group. 

In order to implement the symplectic method, the orig- 
inal second order Lagrangian in the velocity, given in 
(14. ip . is reduced into a first-order form, given by, 



£(0) = n,a^ - V°, 



(4.2) 



with 



2 ° 



^AK-l)-ia;^ (4.3) 



where prime represent spatial derivatives, respectively. 
The symplectic coordinates are S^a'' = (ca, T^a, A) and the 
index (0) indicates the zeroth-iteration. The symplectic 
tensor given by Eq. ()2.3|) is computed in this case as 



/ -Sab 
/(°^ ^\ Sba \ S{x- y). (4.4) 
\ 

This matrix is singular, thus, it has a zero-mode. 



(4.5) 



whose phase space metric is 



(3.41) 



Contracting this zero-mode with the gradient of the sym- 
plectic potential V'^^-', given in Eq. (|4.3p . the following 
constraint is obtained, 



that has a first class constraint, xi, which generates the 
infinitesimal transformations given in (j3.31l) . 



IV. THE 0(N) INVARIANT NONLINEAR 
SIGMA MODEL 



In this subsection, the hidden symmetry present in the 
0{N) nonlinear sigma model will be disclosed enlarging 
the phase space with the introduction of WZ field via dual 
gauge-invariant formalism. We first apply the symplectic 
method to the original second-class model in order to 
show the second-class nature of the model, and also to 
obtain the usual Dirac's brackets among the phase space 
fields. After that, we unveil the hidden gauge symmetry 
of the model which dwells on the original phase space. 

The 0{N) nonlinear sigma model in two dimensions is 
a free field theory for the multiplet aa = (ui, CT2, . . . , (T„) 



(4.6) 



In agreement with the symplectic formalism, this con- 
straint must be introduced into the canonical sector 
of the first-order Lagrangian (|4.2p through a Lagrange 
multiplier p and then, we obtain the first-iteration La- 
grangian as 



with 



(4.7) 



(4.8) 



The symplectic coordinates are ^i^' = (ctq, TTa, p) with the 
following one-form canonical momenta. 
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= 0, (4.9) 
The corresponding symplectic tensor /^^^ given by, 



-Sab o-a 

6ab I 6{x-y), 
-Gb 



(4.10) 



is singular, thus, it has a zero-mode that generates a new 
constraint. 



^1 



(4.11) 



Introducing the constraint ^i into the first-iteration 
Lagrangian (|4.7p through a Lagrange multiplier the 
second-iteration Lagrangian is obtained as 



£(2) = vr.a'^ + p{al - l) + ((^a^'') - 



(2) 



(4.12) 



with V^(2) ^ ^(1) I 

sii=o- The enlarged symplectic co- 

ordinates are fa — {aa,T^a, pX) and the new one- form 
canonical momenta are 



4!^= -a, 

= - 1 
4'^ =^a^^ 



The corresponding matrix /^^^ is 



—Sab CTa TTa 

(^ba CTq 



-TTf, -CT;, 



5{x - y), 



(4.13) 



which is a nonsingular matrix. The inverse of Z'-^-' furnish 
the usual Dirac brackets among the physical fields, given 

by, 

{aa{x),ab{y)Y = {7ra(x), 7rf,(y)}* = 0, 
{aa{x),iTb{y)Y = (Sab - ^) S{x - y), (4.14) 

{Tra{x),TTb{y)}* = (^aTTb J^bT^".) _ 

This means that the NLSM is not a gauge invariant the- 
ory. 



At this stage we are ready to implement our proposal. 
In order to disclose the hidden symmetry present within 
the NLSM, the original phase space will be extended with 
the introduction of WZ field following the symplectic 
gauge-invariant formalism. This process is based on the 
introduction of two arbitrary functions, ^'(a'a, T^a, and 
G{aa,Tra,0), into the first-order Lagrangian as follows. 



£(0) = naa" + "ife - V^°\ 
where the symplectic potential is 



(4.15) 



V^°^ - ^ + 1 A (a^ - 1) + i + G{aa,na, 9), (4.16) 

with G{aa,Tra,0) satisfying the relations given in 
Eqs.dHII) and ((2TT21) . 



The symplectic coordinates are 



(0) 



(fa, TTa, A, ( 



with the following one-form canonical momenta, 

il°)-o, 



A 



(0) 



1 



1), 



= 0. 



(4.17) 



As established by the symplectic gauge-invariant for- 
malism, the corresponding matrix f^'^\ given by 



m 



/ 





-Sab 







Sba 



















V 










9-^1 





Six-y), (4.18) 



must be singular, this fixes the dependence relations 



of arbitrary function 5*, namely. 



0, i.e, * 



\I'(crQ, tTq, 0). This matrix has a zero-mode, identified as 
being the gauge symmetry generator. To pull out the 
hidden symmetry, this zero-mode must satisfy the rela- 
tion p.lSp . allowing then the computation of ^t. 

Let us start considering the symmetry generated by 
the following zero-mode. 




(4.19) 



Since this zero-mode and the symplectic matrix (|4.18p 
satisfy the relation ()3.18p . ^ is determined as 



(4.20) 
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where "c" is a constant parameter. This completes the 
first step of our formahsm. 

The second step begins with the imposition that no 
more constraints are generated by the contraction of the 
zero-mode with the gradient of the potential. The correc- 
tion terms as functions of 9 can be explicitly computed. 
The first-order correction term in 9, t/^^-*, determined af- 
ter an integration process, is 

g<-'\aa,na,9) = -aa7:ae. (4.21) 

Substituting this expression into Eq. (|4.16p . the new 
Lagrangian is 

(4.22) 

However, the model is not yet gauge invariant because 
the contraction of the zero-mode v^^^ with the gradient 
of the potential V'^ produces a non null value, indicating 
that it is necessary to compute the remaining correction 
terms ^/^"^ as functions of 9. It is carried out just de- 
manding that the zero-mode does not generate a new 
constraint. It allows us to determine the second order 
correction term Q^^\ given by 

g^^^ = +l<^lO^- (4.23) 

Substituting this result into the first-order Lagrangian 
(|4^ . we have that, 

+ aa7T''9 - ^ay. (4.24) 

Now the zero-mode j^^^-' does not produce a new con- 
straint, consequently, the model has a symmetry and, in 
accordance with the symplectic point of view, the gener- 
ator of the symmetry is the zero-mode. Due to this, all 
correction terms ^y*-"-* with n > 3 are zero. 

At this moment, we are interested in recovering the in- 
variant second order Lagrangian from its first-order form 
given in Eq. (|4.24l) . To this end, the canonical momenta 
must be eliminated from the Lagrangian (j4.24p . From 
the equation of motion for tTq, the canonical momenta 
are computed as 

TTa = fTa + (Ta9. (4.25) 

Inserting this result into the first-order Lagrangian 
(|4.24p . we have the second order Lagrangian as 

C=^ d^aad^a^ - {daa^)9 l) A, (4.26) 



with the following gauge invariant Hamiltonian, 

'H=\irl + \a'l-{aai^'^)9+\\{al-l) + \ale\ (4.27) 

Both Lagrangian (|4.26p and Hamiltonian (j4.27p are gauge 
invariant. From the Dirac point of view, Jli arises as a 
secondary constraint from the temporal stability imposed 
on the primary constraints, tta and irg, and plays the role 
of the Gauss law, which generates the time independent 
gauge transformation. 

To proceed the quantization, we recognize the states 
of physical interest as those that are annihilated by VLi. 
This gauge invariant formulation of the NLSM was also 
obtained by one of us in [2^ with the introduction of WZ 
fields, as established by the iterative method ITll . and by 
another authors using the BFFT formalism [30|. 

As established by the symplectic formalism, the zero- 
mode is identified as being the generator of the infinites- 
imal gauge transformations 5S,a'^ — ev^^\ namely, 

5(7a = 0, 

5\ = 0, (4.28) 
59 = e. 

For the transformation above the Hamiltonian changes 
as 

5H = 0. (4.29) 

Similar results were also obtained in the literature us- 
ing different methods based on Dirac's constraint frame- 
work [29l - [33 |. However, these techniques are affected 
by some ambiguity problems, as said before, that nat- 
urally arise when the second-class nature of the set of 
constraints turns into first-class with the introduction of 
the WZ fields. In our procedure, this kind of problem 
does not arise and consequently the arbitrariness disap- 
pears. 

Henceforth, we are interested in disclosing the hidden 
symmetry of the NLSM lying on the original phase space 
(cto, TTq). To this end, we use the Dirac method to obtain 
the set of constraints of the gauge invariant NLSM de- 
scribed by the Lagrangian (|4.26l) and Hamiltonian (|4.27p , 
given by, 

h = -\{^l-l), (4.30) 

and 



^2 = (TaT^a ~ CrlO , (4.31) 
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where tt\ and ng are the canonical momenta conjugated 
to A and 6, respectively. The corresponding Dirac ma- 
trix is singular. However, there are nonvanishing Poisson 
brackets among some constraints, indicating that there 
are both second-class and first-class constraints. This 
problem is solved separating the second-class constraints 
from the first-class ones through constraint analysis. The 
set of first-class constraints is 



X2 = -^(a^-l)+^e, (4.32) 
while the set of second-class constraints is 



?2 = (JaTTa - ^a^- (4.33) 

Since the second-class constraints are assumed to be 
equal to zero in a strong way, and using the Maskawa- 
Nakajima theorem (s^l, the Dirac brackets are con- 
structed as 



Saa = e{aa,X2} = 0, 
SiTa = e{7rQ, X2} = ea-a, (4.39) 
6X = 0. 

where e is an infinitesimal parameter. It is easy to ver- 
ify that the Hamiltonian (j4.35p is invariant under these 
transformations because aa are eigenvectors of the phase 
space metric (M^) with null eigenvalues. In this section 
we reproduce the results originally obtained in [36|] from 
an alternative point of view. 

V. THE GAUGE INVARIANT BOSONIZED 
CHIRAL SCHWINGER MODEL 

It has been shown over the last decade that anoma- 
lous gauge theories in two dimensions can be consistently 
and unitarily quantized for both Abelian 0, 113, HI] and 
non-Abelian fSy, '40'| cases. In this scenario, the two di- 
mensional model that has been extensively studied is the 
CSM. We start with the following Lagrangian density of 
the bosonized CSM with a > 1, 



{ai{x),(jj{y)y 



0, 



{ai{x),Trj{y)}* = S^jd{x-y), 



(4.34) 

Hence, the gauge invariant Hamiltonian is rewritten as 



i, 1 2^ 1 ,2 l(aavr«)2 ^ A,_^2 

H = — 7r„ H — (7 „ 

2 " 2 ° 2 (Taa" 

= ^TT,M,jTT^ + ^(t'1 + ^{al - 1), 
where the phase space metric Mij , given by 



^2 ' 



(4.35) 



(4.36) 



which is a singular matrix. The set of first-class con- 
straints becomes 



Xi = ""A, 

X2 - - !)• 



(4.37) 



Note that the constraint X2, originally a second-class con- 
straint, becomes the generator of gauge symmetries, sat- 
isfying the first-class property 



q^aA^Ar 



(5.1) 



Here, F^^, = d^Ai, -d^Afj,, g^i. = diag(+l, -1) and e"^ = 
—e^'-* = eio = 1. The symplectic method will be used here 
to quantize the original second-class model. Then, the 
Dirac brackets and the respective reduced Hamiltonian 
will be determined as well. In order to implement the 
symplectic method, the original Lagrangian written in 
second order in velocity, given in (|5.ip . is reduced into 
its first-order as follows, 



(5.2) 



where the zeroth-iterative symplectic potential U^^^ is 



U^'^ = + r) A, K + \e{a ~ I) A, 



97^0 - 



\q\a + l)A, 



(5.3) 



where dot and prime represent temporal and spatial 
derivatives, respectively. The zeroth-iterative symplectic 
variables are s}a^ = ((/), tt^, Ag, Ai, tti) with the following 
one- form canonical momenta A^, 



{X2,H} = Q. 



(4.38) 



Due to this result, the infinitesimal gauge transforma- 
tions are computed as 



(0) 
(0) 



(5.4) 



4"^ = Af = 4") = 0. 

T^d, A.Q TTi 
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the zeroth-iterative symplectic tensor is obtained as 



/O -1 \ 
1 

0-1 

Vo 1 / 



5{x - y). (5.5) 



This matrix is obviously singular, thus, it has a zero- 
mode that generates a constraint when contracted with 
the gradient of the potential C/*^°\ given by. 



AO) 



Tr[ + q'^ia - l)Ao + q'^Ai + q-n^ + qcf)' , (5.6) 



that is identified as being the Gauss law, which satisfies 
the following Poisson algebra. 



(5.7) 



Substituting the constraint J7i into the canonical sec- 
tor of the first-order Lagrangian through a Lagrange 
multiplier rj, we have the first-iterative Lagrangian L'<^\ 
namely, 



i« =7r^(^ + 7riii+Oi?7-C/W, (5.8) 
with the first-order symplectic potential given by 



+ q^Ax + qiT4, + q(j)'] 

- Ai[-q'K^-]^q^{a + l)Ai 



(5.9) 



where U'^^^ = U^'^\ Therefore, the symplectic variables 
become = {(t),'Kfj,,Ao,Ai,TTi,rj) with the following 
one-form canonical momenta. 



.(1) 
Ml 



A\: =7ri. 



A«=A«=4\)=0, (5.10) 
The corresponding matrix f'^^) is then 





1 





-qdx 



-1 







-q 









-q\a-l) 









qdy 


\ 








q 








q\a-l) 







-1 






1 





dy 




-q' 


-dx 





/ 



which is a nonsingular matrix. The inverse of f^^\x,y) 
furnishes, after a straightforward calculation, the Dirac 
brackets among the physical fields. 



r 



5{x - y), 



(5.11) 



{^{x),<Piy)y 
{(l){x), 7r^(y)}* 

{ct>{x),Ao{y)r 

{ct>{x),A,{y)r 
{(l){x) , 7ri(y) }* 
{'K4,{x),n^{y)}* 

{'K^{x),Ao{y)}* 

{'K^{x),Ai{y)}* 

{'K^{x),'Ki{y)}* 

{A,{x),Ao{y)r 

{Ai{x),A^{y)r 
{A,{x),My)r 

U,{x),Ao{y)r 
{7ri(a;),7ri(y)}* 



0, 

5{x-y). 
1 



q{a-l) 



5{x-y), 



0, 
0, 
0, 



1 



0, 
0, 

0, 

5{x - y) , 
0. 
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This means that the model is not a gauge invariant the- 
ory. 

The gauge symmetry inside the model will be disclosed 
via a new gauge- invariant formalism that does not require 
more than one WZ field. The fundamental concept be- 
hind the symplectic gauge-invariant formalism dwells in 
the extension of the original phase space with the intro- 
duction of two arbitrary function ^'(0, tt^, Aq, Ai, tti, 0) 
and G((/(), 770, j4o, j4i, TTi, 0), depending on both the orig- 
inal phase space variables and the WZ variable 6, into 



the first-order Lagrangian, right on the kinetic and sym- 
plectic potential sector, respectively. In this way, the 
first-order Lagrangian that governs the dynamics of the 
bosonized CSM, given in Eq. (|5.2p . is rewritten as 

L^"'> ^ Tr^4> + TTiAi + d'f - U^°\ (5.13) 

where 



- Ai[^qn^^ ^q^{a + l)Ai- qcj)'] + G((/), tt^, Aq, Ai, tti, 6*). (5.14) 



The gauge- invariant formulation encompasses two The corresponding symplectic matrix f^^^ reads 
steps: one is dedicated to the computation of ^ while 
the other is addressed to the calculation of G. 

The enlarged symplectic variables are now ^i"'' = 
(0, TTcf,, Aq, Ai,TTi, 9) with the following one-form canoni- 
cal momenta 

= = = 0, (5.15) 



?(o) 



( 





-1 













1 













































-1 













1 



















d(t>y 


9-^1 









(5(x - y). 



(5.16) 



where 



(5.17) 



with EE 6'(a;), Oy = 6{y), = *(a;) and ^-j, = ^'(y). 
Note that this matrix is singular since §^ = 0. Due to 
this, we conclude that ^' = ^'(0, tt^, Ai, tti, 0). 

To unveil the gauge symmetry hidden inside the model, 
we assume that this singular matrix has a zero-mode 
[v'^^'^) that satisfies the following relation. 



^'£\x)f^l{x-y)dy^Q. (5.18) 



From this relation a set of equations will be obtained 
and consequently, the arbitrary function 4' can be deter- 
mined. We will now investigate the symmetry related to 
the following zero-mode. 



{q -qda: 1 -q^ -l) 



(5.19) 
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with bar representing a transpose matrix. 

To start, we multiply the zero-mode (|5.19l) by the sym- 
plectic matrix ()5.16|) . as shown in equation (|5.18|) . Due 
to this, some equations arise and after an integration ^ 
is determined as 



by the arbitrary function G, given in Eq. (|2.1ip . will 
be computed. It is achieved just imposing that no more 
constraints arise from the contraction of the zero-mode, 
given in Eq. (|5.19p . with the gradient of the symplectic 
potential, 



+ qc/)' + qn^ + q^Ai. 



(5.20) 



The corresponding symplectic matrix (j5.16p is rewritten 

as 






-1 











qdy \ 


1 














q 
































-1 













1 





dy 


qdx 


-q 





-q^ 


-dx 


) 



S{x~y) (5.21) 



which is obviously singular. Consequently, it has a zero- 
mode that, by construction, is given by equation (|5.19l) . 

Now we start the second step of the method reformu- 
lating the model as a gauge invariant theory. At this 
stage, the correction terms as functions of 0, embraced 



(0) 



(5.22) 



The first-order correction term in 9, Q^^^ , is determined 

by, 

(0,7^0, ^1,7^1,^0,0) = ^ni0 + q^ia-l)A[0-q^0ni, 

(5.23) 

after an integration. Substituting this expression into the 
equation (j5.13l) . the new Lagrangian is obtained as 

L(°) = 7r00 + TTiii + - C/(°\ (5.24) 

with 



= 1(^1+^1 + 0") - ^0 [< + ^q^ia - l)Ao + q^Ai + qir^ + #'] 



(5.25) 



The Lagrangian in (|5.24p is not yet gauge invariant be- 
cause the zero-mode j/^^' still generates new constraints, 
given by 



/I) 



(0) 



q^{a-l)9" -q^{a-l)9 + q^9, (5.26) 



indicating that it is necessary to obtain the rest of the cor- 
rection terms t/*^"' as functions of 9. This is accomplished 
just imposing that no more constraints are generated by 
the contraction of the zero-mode with the gradient of ex- 
tended symplectic potential. It allows us to determine 
the second order correction term Q^^^ as. 



Hence, the first-order Lagrangian (|5.24l) becomes 

L^") =7r^(/) + 7riii+^'^-J7(°), (5.28) 
with the new symplectic potential 



-q^a-l)9 + q^{a-l)9" + q^ 

a 1 



89 



0, 



S''^--l q'ia~l)9" + lq'9'-lq\a^l)9\ 

(5.27) 
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^Ai[-qT:^- + i)Ai - q(b'] - ^^9 + q^{a ~ 1)9 A[ ~ q^TTi 

qHa- 1)9'' + \q'9^ - \q\a l)9\ (5.29) 

I 



The contraction of the zero-mode v'^'^^ with the new sym- 
plectic potential above does not produce a new con- 
straint. Consequently, the model has a symmetry and 
this zero-mode is the generator of the infinitesimal gauge 
transformations. Due to this, all correction terms C/*^"^ 
with n > 3 are zero. The infinitesimal gauge transforma- 
tions generated by the zero-mode {5£,i = ev^^^) are 

(50 — qe, 
Sn^ = qe\ 
Mo = e. 

Ml = -e', (5.30) 
S-Ki = — g^e, 
59 = -e. 

It is easy to verify that the Hamiltonian, identified as be- 
ing the new symplectic potential U^^\ is invariant under 
these infinitesimal gauge transformation above, namely. 



The time stability condition for the constraint xi just 
allows us to determine the Lagrange multiplier A3 . In this 
way, the gauge invariant version of the model has three 
constraints, namely, ipi, (p2 and xi- The corresponding 
Dirac matrix, given by, 

/ -q'{a~l) \ 
C{x-y)^ \ q^a-l) q^{a-l) \S{x-y), 
\ ~q^{a^l) / 

(5.36) 

is singular. As the Dirac matrix is singular, the model 
has both first and second-class constraints. Through a 
constraint combination, we obtain a set of first-class con- 
straints such as, 

XI = -TTe + * - TTo, (5.37) 
and a set of second-class constraints, given by 



(5J7(°)=0. (5.31) 

At this point, we are interested in the analysis of this 
result and we would like also to demonstrate that the 
anomaly was canceled. It will be carried out from Dirac's 
point of view. From the Lagrangian (|5.28p the chains of 
primary constraints are computed as, 

Xi = -^e + *- (5.32) 

These primary constraints will be introduced into the 
Hamiltonian through Lagrange multipliers. In this way, 
the primary Hamiltonian, 

US-Lary^U^°^ +^IVI+^IX1- (5-33) 

Since the constraint ipi has no time evolution, the follow- 
ing secondary constraint is required 

^2=^i-q\a-l)9, (5.34) 

and no more constraints arise from the temporal stability 
condition. In this way, the total Hamiltonian is 



ip2^^i-q'{a~l)9. (5.38) 

It is again easy to verify that Xi is a first-class constraint, 
identified as the Gauss law, while the others are second- 
class constraints. Note that the anomaly was removed. 
Hence, the Gauss law is also recognized as being the gen- 
erator of the gauge transformation given in Eq. (|5.30l) . 

At this stage, we will compute the degrees of free- 
dom of the gauge invariant model proposed by us. The 
model has one first-class and two second-class constraints 
and the phase space dimensions sum eight dependent 
fields, i.e., (0, tt^, Ai, tti, Aq, ttq, 6*, tt^). The first-class 
constraint eliminates two fields, while the second-class 
constraints eliminate two fields, summing then four fields 
eliminated. Hence, the model has four independent fields, 
i.e., there are two independent degrees of freedom. 

In order to obtain the Dirac brackets, the set of second- 
class constraints, ipi and (p2 , will be assumed equal to 
zero in a strong way. After a straightforward computa- 
tion, the Dirac brackets among the phase space fields are 
obtained as 



U^li = U'-"^ + AiV'i + X2^2 + 71X1- 



(5.35) 
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{</>(x),0(y)}* 
{(t>{x),n^{y)}* 

{ct>{x),Ao{y)r 

{<t>ix),A,iy)r 
{</>(x),7ri(y)}* 

{m,0iy)r 

{TT4x),Ao{y)}* 

{TT^{x),Tr^{y)}* 
{Tr^ix),Ai{y)}* 
{Tr^{x),Tri{y)}* 

{TT^ix), 9{y)r 
{TT^{x),Trg{y)}* 

{A,ix),Ao{y)r 

{A,{x),A,iy)r 
{A,{x),My)r 

{A,{x),9{y)r 

{Ai{x),7Tg{y)}* 

{7r,{x),Ao{y)r 

{7ri(a;),7ri(y)}* 

{Mx),9{y)r 
{TTi{x),TTe{y)}* 
{e{x),Ao{y)r 

{e{x),My)r 

{TTe{x),Ao{y)}* 
{TTg{x),TTe{y)}* 



0, 

S{x - y), 
1 



qia-1) 



S{x - y), 



0, 
0, 
0, 
0, 



1 



0, 
0, 
0, 
0, 
0, 



-dyS{x - y), 



(5.39) 



1 



q^a-l) " 

0, 

5{x ~ y), 
0, 
0, 

1 



dyS{x - y), 



r 



(« 
0, 
0, 
0, 
0, 

5{x - y), 
5{x ~ y), 
0. 



S{x - y), 



Note that the Dirac brackets among the original phase 
space fields were obtained before. After this process, the 
model now have only one first-class constraint, 



X = Xilvi=V2=o = -TTe + * 



(5.40) 



identified as the Gauss law, that satisfies the Poisson al- 
gebra. 



(5.41) 



In this way, the anomaly was eliminated, the symmetry 
is preserved, and the fundamental brackets among the 
original phase space fields were reobtained. Note that 
the Gauss law is the generator of the gauge symmetry 
given in (|5.30p . 

Once more, the number of the independent degrees of 
freedom matches with the result obtained in the second- 
class case. The invariant model has a phase space 



((/), TT^, Ai, TTi, 0, TTg), totaling six dependent fields, and 
has a first-class constraint which eliminates two fields. 
Consequently, the model has two independent degrees of 
freedom. 

At this point, we are interested in commenting about 
the consistency of the gauge invariant version of the 
bosonized GSM. To do this, the remaining symmetry will 
be eliminated with the introduction of the unitary gauge- 
fixing term, given by. 



e = o. 



(5.42) 



Hence, both the noninvariant Hamiltonian and the cor- 
responding Dirac brackets computed in the beginning of 
this section are reobtained, then recovering the anomaly. 
In this way, we conclude that the new symplectic gauge- 
invariant formalism does not change the physics contents 
present inside the model. 



VI. DUALITY FOR THE NON-ABELIAN 
EXTENSION OF THE PROCA MODEL 

The non-Abelian extension of the Proca model has its 
dynamics governed by the following Lagrangian density. 



with 



(6.1) 



(6.2) 



where the antisymmetric tensor C^^ (C^^ = — C"^), are 
in fact a set of real constants, known as the structure 
constants of the gauge group, and satisfy the following 
property. 



Ca r^d I r^a r^d i r^a r^d r\ 
bc^ae "I" '^fh'^nr "T '-'rp'-^nh — 



^eb ac 



ce ab 



(6.3) 



Since we are interested in analyzing the non-Abelian 
Proca model from the symplectic point of view, the La- 
grangian will be reduced to its first-order form as follows. 



-\m'AlAl~\F^^F^^, 

where 

^a^d.^^-gCl^^lA^^m'Al. 
The symplectic variables are given by 



(6.4) 



(6.5) 
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agreement with the symplectic method, this constraint is 
introduced into the kinetic sector of the first-order La- 
grangian through a Lagrange muhipher, namely, 



and the symplectic matrix is 




I (5(3'(^-y)- 

0; 



(6.6) 



Since this matrix is singular, it has a zero-mode that 
generates the constraint ila, given by equation (|6.5p . In 



J 



(6.7) 



The new group of symplectic variables is — 
{Af , Trf , Aq , f]'^) , and the new symplectic matrix is 





^ba—d ( 



\gC'/7r^ix) 



-S,.d'- -gCS'Trfiy) \ 

S'^^'d^ - gCfAfiy) 

m^S"'' 

-S^^d"!; + gCfAj{x) -m^5^'' / 



(6.8) 



This matrix is nonsingular and its inverse leads to the 
commutation relations among the dynamical variables, 
given by 



{A1{x),A]{y)}^Q, 
{A1{x),^]{y)}^5'^H,,5^^\x-y), 
K(x),^^^(y)} = 0, 

{A1{x),Al{y)} = -1^5'^^d^5'^^\x-~y) 

-l,^CfA'l{x)5{x--y), 

{AUx),A',iy)} = -2gm'CfAl{x)5^^^\x~y) 
{K{x),Al{y)} = -4cf <(x)<5(3)(f -y). 



(6.9) 



vI/, = vI/,(A,^7rf,Ag,r), 

OO 

Gee G(A^<,Ag,0°) = ^ ^"(A^, <, AJJ, ^ ) 



n=0 



(6.11) 



and where the G function obeys a boundary condition 
given by. 



G EE (A^<, = 0) = g"(A^<, = 0) = 0. 

(6.12) 

In this context, the corresponding symplectic matrix is 



This completes the analysis of the noninvariant descrip- 
tion of the model. 

The model now can be reformulated as a gauge in- 
variant field theory. This will be carried out in the 
context of the symplectic gauge-invariant formulation. 
In agreement with this formalism, the first-order La- 
grangian (j6.4p can be rewritten as 



- \m'A1^^ - \m^AlAl - -F^.F^, - G, 



(6.10) 



where the arbitrary functions are 



/ 



/ 



(0) 







S^'Hx-y). 



dASix) 

\ dA^^iy) d7T''{y) aAg(y) 

(6.13) 

In order to determine the functions, we analyze the 
symmetry related to the following zero-mode. 



p{0) ^ (^gx 1) 



(6.14) 



with df = which satisfies the following condition, 



i>^''\x)fo.p{x-w)^Q. 



(6.15) 



This condition produces a set of differential equations 
which allows us to compute the function as 
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= -diiriix). (6.16) 
Consequently, the first-order Lagrangian is rewritten as 



wliere the symplectic potential is 



(6.17) 



(6.18) 



It completes the first step of the symplectic gauge- 
invariant formulation. 

To unveil the hidden symmetry inside the model, the 
zero-mode i^^^^ does not generate new constraints, conse- 
quently, we have the following relation. 



(6.19) 



From this relation we can compute the whole set of cor- 
rection terms as functions of 9°'. The linear correction 
term in 6 is computed as being. 



dAi{w) d9f{w) 
After an integration we have that 



0. 



(6.20) 



-gCldf{A-,{x)4ix))efix) - m''{dfA'f)ef{x) 



(y). (6.21) 

Now, we will compute the quadratic term, namely, 
I 



g^'\x) = \f [d^[ /a 



dg'-^Hx) g^Hx)_ 

dA{{w) ^ oefiw) 



0. 



3.22) 



Integrating this relation in 6^{w), the quadratic correc- 
tion term is obtained as 



g^'\x) = l-m'{di9f{x)r 



+ 



\f [ or f [{d!A^j,)9' 

^ Jef{x) Jw Jy 



where 



Ail 



(6.23) 



dA{{w) dA\{y) + ^^-^"^ 



dA{{w)dA\{y)' 

(6.24) 



In this way, two correction terms as functions of 9a 
{g^^\x) and g'^^\x)) remain to be computed. Let us 
compute the first one. It can be done from the following 
relation, 



g^'\x) 



d9g{z) 
and we can write that, 

g^'\x) 



(6.25) 



Finally, the last correction term is. 



0f(x) Jw 



dA^^iv)dAl^iz) 



iv)- 



(6.26) 



Therefore, the gauge invariant first-order Lagrangian is where the gauge invariant Hamiltonian, identified as be- 
ing the symplectic potential, is given by 



£(0) = niA^ - (5i<)r - , (6.27) 
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lo,M I or f [idf^l,)e''{y)] -\[ ^Jdnl dr [ 9f^^/iy) 

^ Jw Jy ^Jes{z)Jz Jef{x)Jw Jy CIAg{Z} 

^ ^ I I dll I or I dy j'^', 0\y). (6.28) 



This completes our proposal. 

At this stage, we would like to disclose the hidden sym- 
metry present inside the model using the Dirac point of 
view. To this end, we start with the set of primary con- 
straints. 



C 



with 



2 red 



(6.33) 



(6.29) 



For the first set of constraints, the time stability con- 
dition is satisfied (f7° — 0). For the second one, the 
following secondary constraints are required. 



Due to this, the total Hamiltonian is 



(6.30) 



(6.31) 



where A;^, C,\ and are Lagrange multipliers. Since the 
Poisson brackets among those constraints are 

{17?(a:),x?(y)} = 0, 

{ni{x),xl{y)} = Q, (6.32) 

{X'i{x),xl(y)} = 2gCfxi{x)6^'\S~y) 

-2gm^CfAt{x)5'^^\x-y), 

no more constraints arise. Notice that some brack- 
ets above are null, indicating that there are both first 
and second-class constraints. Indeed, the first-class con- 
straint is ^1 and the second-class are Xi and X2- 1^ 
agreement with Dirac's procedure, the second-class con- 
straints can be taken equal to zero in a strong way, this 
allows us to compute the primary Dirac brackets. Due to 
the Maskawa-Nakajima theorem [s^, the primary Dirac 
brackets among the phase space fields are canonical. To 
demonstrate this, the brackets are computed explicitly. 
The Dirac matrix is 



The inverse of the Dirac matrix is 



-(5''^ 



3.34) 



(6.35) 



In accordance with the Dirac process, the Dirac brackets 
among the phase space fields are obtained as 



{A'^{x),A]{y)Y 
{A'iixli^^iy)}* 

{A1{xlAl{y)Y 



{nt{x)^Al{y)Y 

{Kix),n'^,{y)r ^0, 

MS(a;),Ag(y)r = - 
1 

{Anx),9''{y)r =0, 

{K{x),e''{y)r = o, 

{A^,{x),9\y)r^0, 
0, 
0, 



0, 

S^''S^^\x-y), 
-^9fS^'\x-y) 

J^CfA{{x)6(^Hx-yh 

-^gCfntS^'Hx-y), 



;CfA^oix)S^^\x-y}, 



(6.36) 



K"(x),4(y)}* 
K(x),4(y)}* 

{AS(x),4(y)}* 

{0%x),4{y)r 
{e'^{x),d\y)r 
W(x),4(2/)}* 



,Cfdf7T9x)S^^\x-y), 



6''''6'^^\x 

0, 
0. 



y), 



Finally, the infinitesimal gauge transformations are ob- 
tained, namely. 
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SA1 = -d^e", 
5n^ = 0, 
SA^, = 0, 
60" = e^ 
= 0, 



3.37) 



which lead to the invariant Hamiltonian. 

To demonstrate that the gauge invariant formulation of 
the non-Abelian Proca model is dynamically equivalent 
to the original noninvariant model, the symmetry is fixed 
by using the unitary gauge fixing procedure, 



= ff" « 0, 
which leads to the bracket below, 

{nnx),v^\y)} = -s^'s^'Hx-y). 



3.38) 



3.39) 



Due to this, a new Dirac brackets must be computed. 
The corresponding Dirac matrix for this set of constraints 
is 



C 



-1 

1 



(6.40) 



Using the inverse of this matrix, the Dirac brackets 
among the physical phase space fields are computed, 
which is equal to the one calculated in the original de- 
scription given by equation (16. 9p . This result demon- 
strate that the symplectic formalism can be seen in fact 
as a dual mapping between the original theory and the 
final one. 



VII. HIDDEN SYMMETRY IN THE 
ROTATIONAL FLUID MODEL 

Recently, some of us have proposed a Wess-Zumino 
(WZ) gauge-invariant version for the isentropic irrota- 
tional fluid model [ZTl. In that work, we have demon- 
strated that the irrotational fluid model has a set of dy- 
namically equivalent WZ gauge invariant versions. Fur- 
ther, in that paper the extra global symmetries, namely, 
Galileo antiboost and time rescaling, first obtained in 
(43 [. were lifted to local symmetries. 

In this section, we propose an investigation of the sym- 
metries of the rotational fiuid model, but now with an 
extra term, like kp [diO + adif3)'^ , where fc is a constant, 
p = p{t,r) is the mass density and 9 = 9{t,r) is the ve- 
locity potential. This term introduces a dissipative force 
into the model. The main motivation is that in the real 
world there is always dissipation, albeit sometimes ex- 
tremely small. Our purpose is, to verify if the extra sym- 
metries found in [41| are broken with the introduction of 



this new term in the model. We will introduce a final 
result that is new in the literature. 

As was demonstrated in [4lj , the isentropic irrotational 
fluid model has a dynamically equivalent family of WZ 
gauge-invariant descriptions. Thus, in order to establish 
our ideas, we are going to investigate how the inclusion 
of a dissipative term affects the dynamics of the fluid 
model considering only one of these zero modes presented 
in [4l|. Note that in the prescription of the symplec- 
tic formalism, for each zero-mode chosen we will have a 
gauge-invariant version for the model. 

In the next subsection, the scalar rotational fluid the- 
ory will be analyzed from the symplectic point of view 
|23| , and the fundamental Dirac brackets of the fields will 
be computed. 



A. Symplectic analysis of the rotational fluid 

In this section, we will analyze the rotational fluid dy- 
namic model from the dual embedding point of view. 

It is well known that systems that have vorticity 
and/or viscosity present Casimir invariants which ob- 
struct the construction of a canonical formalism for fiuid, 
as demonstrated in [43[. However, this obstruction can 
be eliminated using the Clebsch parameters, as shown 
by Lin and by two of us in [45|. In fact, with the 



introduction of the Clebsch parameters, it is possible to 
obtain a Lagrangian density for the rotational fluid with 
dissipation, in 3-dimensional, as being 

C^-pie + aP)~V, (7.1) 

where the symplectic potential is 

y = i(l - k)p{d,e + ad,(3){d'd + ad'13) + V{p). (7.2) 

The symplectic coordinates are = {p, 9, a, (3) with the 
corresponding zeroth-iterative one-form canonical mo- 
menta given by 



= -P. 
^ 0, 
Af = -ap. 

The zeroth-iteration symplectic matrix, given by 



/ 



V 





(5(f-r') 

a 



—S{f— f") —a 




(7.3) 
(7.4) 

(7.5) 



is a nonsingular matrix and, consequently, the model is 
not a gauge-invariant field theory. As settle by the sym- 
plectic formalism [2^, the Dirac brackets of the phase 
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space fields are acquired from tlie inverse of the symplec- 
tic matrix, and are given by 



Mr),0(r')}*=<5(r-r'), 
P 

{a(f),/3(r')}* = i(5(f-r'), 
P 



(7.6) 



while the remaining brackets are null. This completes 
the noninvariant analysis. 

In the next subsection, we will apply the symplectic 
embedding formalism j28| in the 3-dimensional fluid dy- 
namical model with dissipation and, as a consequence, 
we will obtain a gauge-invariant version for this model. 
Although the symplectic formalism does not restrain the 
dimension of the model, we choose a 3-dimensional de- 
scription for the rotational fluid in order to put our work 
in a correct perspective with others. 



B. Obtaining the gauge-invariant version of the 
model 



In this section, we will obtain the WZ gauge-invariant 
version of the fluid theory. We will consider initially a 
general interaction potential and after that, a specific po- 
tential {V = g/ p) will be chosen. Following the prescrip- 
tion of the formalism, two arbitrary functions, ^' and G, 
depending on the original phase space fields and the WZ 
field [rf), must be added to the model. The former is 
introduced into the kinetic sector and the latter one into 
the potential sector of the first-order Lagrangian. The 
process starts with the computation of ^ and finishes 
with the computation of G. 

In order to reformulate the model as a gauge-invariant 
field theory, let us start with the first-order Lagrangian 
C^°\ equation (|7.ip . with additional arbitrary terms 
(*,G), 



(7.7) 



where 



- i(i - k)p{d^e + ad,p){d'e + ad'p) + v{p) + g, 

(7.8) 

and where ^ = ^(p, 9) and G = G(p, 0, 77) are arbitrary 
functions to be determined. Now, the symplectic coordi- 
nates are S^^^ = (p, 9, a, /3, 77) while the symplectic matrix 
is 






-5{f-r') 





—a 


5[f-r') 




















-p 


a 





P 







5*P 






Sp(r') 


5e(r') 


5a{r' ) 


5/3 (r' 



5*-, 



5a (r) 

WW) 




where '^p= '^{p{r),9{r)) and -^f, = ■^{p(r'),9{r')). 



As established by the dual embedding method, the cor- 
responding zero-mode i'^"^ (r) satisfies the following rela- 
tion 



dfi>(0)«(f) /,-^(f,r') = 0, 



(7.9) 



which produces a set of equations that allows to deter- 
mine ^ . At this point, it is very important to notice that 
the formalism unveils the U{1) hidden gauge symmetry 
of the physical model because the zero-mode does not 
generate a new constraint. Considering a general zero- 
mode, 



^ ( a h c d ~l 



(7.10) 



we have the following set of differential equations 
j df (^5{f- + da5{f- r')^ = 0, 

After a direct calculation, we obtain 

51(7") = -fop - adp + a9 + aa/3 + cpP. 



(7.11) 



(7.12) 



In order to have a solution, we consider that a — c — 
then 



and 



,>(o) = (0101-1 



*(f) = -{l + a)p. 



(7.13) 



(7.14) 



Now, we begin with the second step of the method to 
reformulate the model as a WZ gauge-invariant model. 
The zero-mode j>(°' does not produce a constraint when 
contracted with the gradient of the symplectic potential, 
then. 



df {(1 + a)p{d[9 + ad[l3)d[6{f - f) 



^ \ d9 

71 — 1 ^ 



)} 



= . (7.15) 



This expression produces a general differential equation, 
that allows the computation of all the correction terms 
as functions of 77 enclosed into G(p, 0, rj). To compute the 
first correction term as function of 77, Q^^^ , we pick up the 
terms in equation (j7.15l) with zeroth-order in rj, thus 



df \{l + a)p{d'i9 + adl(3)d'i5{f -f)- 



drj 



= 



(7.16) 
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where = After a straightforward calculation, the write 
linear correction term as function of rj is obtained as 

= (1 + a)p (9,0 + a9,/3) a,77. (7.17) 

For the quadratic correction terms in equation (I7.15|) . we 
have that 



ag^i) 9^(1) 



de 



dfi 



drj 



(7.18) 



After a direct calculation, the second-order correction 
term is obtained as 



g(2) = (1+^)' 



pdiTjd'-ri. 



(7.19) 



9^(2) ag(2) ag(3) 



de 



dfi 



dr] 
dr] 



0, 
0, 



(7.20) 



which allows to conclude that ty*^"^ = for n > 3. Hence, 
the gauge-invariant first-order Lagrangian is written as 

= -p{0 + a/3) - (1 + a)pr) - V^°\, 
= -p(^ + 77) - ap{l3 + //) - V^(°) 

where the symplectic potential is 



(7.21) 



For the cubic correction terms in equation (j7.15l) , we can 



= 77(1 - k)p{d,9 + ad,J3){d'9 + ad' (3) + (1 + a)p (9,;i 



ad^p) d,rj + itt^pd^T^d'rj + Vip), 



(1 - k)p [d,{e + r/) + ad,{P + r/)] [d'{9 + r?) + ad'{P + tj)] + V{p) . 



(7.22) 



To complete the gauge-invariant reformulation of the 
model, we will compute the infinitesimal gauge transfor- 
mation. In agreement with the method, the zero-mode 
is the generator of infinitesimal gauge transforma- 
tions {SO = £!>("'). Then, 

6pir,t)^0, 

S9{r,t) = e{f,t) Sif-f'), 

da{f,t) = 0, (7.23) 

5/3ir,t) =e{r,t) 5{f-r'), 

5'q{r,t) = —s{r,t) S{r~f'), 

where e(r, t) is an infinitesimal time-dependent param- 
eter. In fact, under the infinitesimal transformations 
above, the invariant Hamiltonian (V"*^")) changes as 

(5f(°) = 0. (7.24) 
Considering the following transformations, 

P^P + ri, (7.25) 

then the Lagrangian density, equation (|7.2ip . and the 
Hamiltonian, equation (|7.22p . become 

= -p(^ + a^) - (7.26) 
V^^) = ]^p{di9 + adrP){d'9 + ad'P) + V{p) . 

These expressions are similar the original expressions for 
the Lagrangian in equation (j7.1l) . and the Hamiltonian 
in equation (|7.2p . respectively. 



Thus, at this point, it is important to point out that 
exist a hidden symmetry into the rotational fluid model. 
Note that, in Section [TTl it was not possible to realize this 
symmetry. Another interesting feature discovered here is 
that for the set of the differential equations obtained we 
have no other solutions, i.e., the model has, in fact, only 
one hidden symmetry. Based on the investigation done 
by some of us in Ref. (4l| , where the extra global symme- 
tries proposed in are lifted to the local status, we can 
conclude that these extra global symmetries do not exist 
in the rotational fluid model. Some other considerations 
will be depicted in the next section. 

VIII. FINAL DISCUSSIONS 

From the Dirac point of view, a system classified as 
a gauge invariant theory is one that has only first-class 
constraints. When a theory has second-class constraints, 
the gauge invariance can be recovered by converting the 
second-class constraints into first-class constraints. In 
the literature there is a great variety of techniques, with 
pros and cons, to promote this kind of conversion. 

In this work we are concerned not only with the gauge 
invariance of second-class systems but also in obtaining 
a theory dual equivalent to the original one. We be- 
lieve that dual embedding formalism is the most ade- 
quate technique because it is not affected by ambiguity 
problems related to the introduction of the WZ variables 
[13| . Besides, as demonstrated by some of us, this method 
has the advantage that a convenient choice of a conve- 
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nient zero-mode can lead to a theory dual equivalent to 
the parent one through the elimination of the WZ terms 
[J, [2|. This means a new interpretation of the method. 
However, as recently demonstrated by some of us |47i| 
the choice of the zero-mode must obey some "boundary 
conditions". In other words, we can say that the phys- 
ical coherence must guide us to choose the correct (or 
convenient) zero-mode. But it is not an impediment to 
obtain, as mentioned just above, a whole family of dual 
equivalent actions. 

Speaking in another way, we can say that this particu- 
lar mapping between the parent action and the respective 
final gauge invariant theory can be interpreted as a kind 
of duality. The gauge invariance of the final actions ob- 
tained here was demonstrated via the Dirac analysis. The 
characterization of the final action as a first-class system 
corroborates the success of this process of "dualization" . 

Firstly in this paper we used a kind of toy model, the 
Proca model, to illustrate the procedure. Although be- 
ing a toy model, the resulting "dual" action is new. Af- 
ter that, we apply the formalism to the non-linear sigma 
model and to the chiral Schwinger model. In the NLSM, 
a hidden symmetry lying on the original phase space was 
disclosed, oppositely to other approaches [H, IH, [30| . 
where the symmetry resides on the extended WZ phase 
space. 

In the CSM, the chiral anomaly was eliminated and the 
gauge symmetry was recovered. It is important to notice 
that this result was achieved introducing one WZ field 
while other schemes thrive with the introduction of two 
or more WZ fields, which is the origin of the ambiguity 
problem. 

Besides, we showed in the context of a non-Abelian 
model (the non-Abelian Proca model) that the dual em- 
bedding formalism can be used without any restrictions 
with the algebra obeyed by the noninvariant model, while 
other constraint conversion techniques work since the al- 
gebra was previously and necessarily taken into account. 

We also have proposed a gauge-invariant version for 



the rotational fluid model. As a consequence, we have 
demonstrated that the hidden symmetry found is unique. 
Although we have studied the rotational fluid model 
with an extra term, which introduces dissipation into 
the model, the results are also valid without viscosity 
{k = 0). We have noted that, although having dissipa- 
tion, this fluid model presents hidden symmetry, which 
does not belong to the other symmetries group obtained 
for the irrotational fluid model [4l|. Due to this, the lo- 
cal version of the extra global symmetries does not 
exist in the rotational fluid model, with dissipation or 
not. Furthermore, the physical meaning of the hidden 
symmetry can be interpreted. Consider a flow of fluid 
with viscosity in a tube. It is well known that there 
are many layers j46| . each layer flowing with a specific 
velocity: the velocity of the layer at the center of the 
tube is maximum and, distant from the center, the ve- 
locity decreases in accordance with the increasing of the 
distance from the center. It also happens when vortices 
are present. The are many layers in a vortex and the 
one that is near to the center of the vortex has increased 
its velocity. Due to this, the velocity of each layer can 
be transformed into the velocity of the other layer by 
using the equation s in equation (|7.25|) . Therefore, the 
dynamics of the fluid is preserved and governed by both 
Lagrangian and Hamiltonian which present the hidden 
symmetry. 
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